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ABSTRACT 



We present simultaneous spectropolarimetric observations of four visible (630 nm) and three infrared (1565 nm) spectral lines from 
the German Vacuum Tower Telescope, together with speckle-reconstructed filtergrams in the G band and the Ca n H line core from 
the Dutch Open Telescope. After alignment of the data sets, we used the G-band intensity to locate bright points (BPs) in the moat 
of a regular sunspot. With the cospatial and cotemporal information provided by the polarimetric data, we characterize the magnetic, 
kinematic, and thermal properties of the BPs. We find that (a) 94% of the BPs are associated with magnetic fields; (b) their field 
strengths range between 500 and 1400 G, with a rather flat distribution; (c) the contrast of BPs in the G band depends on the angle 
between the vector magnetic field and the line of sight; (d) the BPs harbor downflows of magnetized plasma and exhibit Stokes V pro- 
files with large area and amplitude asymmetries; (e) the magnetic interior of BPs is hotter than the immediate field-free surroundings 
by about 1000 K at equal optical depth; and (f) the mean effective diameter of BPs in our data set is 150 km, with very few BPs larger 
than 300 km. Most of these properties can be explained by the classical magnetic flux tube model. However, the wide range of BP 
parameters found in this study indicates that not all G-band BPs are identical to stable long-lived flux tubes or sheets of kG strength. 

Key words. Sun: magnetic fields, spectropolarimetry 



1. Introduction 

When observed at sufficient angular resolution, the solar pho- 
tosphere shows small elements that appear bright in continuum 
intensity. Their contrast is particularly high in the G band, the 
wavelength region around 430 nm dominated by absorption lines 
of the temperature-sensitive CH molecule, and other molecu- 
lar bands such as the CN band head at 388 nm (Zak harov et alj 
l2005HUitenbroek & Tritschlerfl2006l) . The discovery o f these G- 
band bright points (BPs) dates back to lMulletj ([1983), although 
they can be iden tified with the " bright poin ts" or "filigree" de - 
scribed earlier bv lDunn & Zirketj d 19731) and lMehltreltetl d 19741) . 

The BPs are observed preferentially along intergran- 
ular lanes, both in active regions and in the quiet Sun. 
Today it is believed that they trace magnetic flux concen- 
trations. The relationship between BPs and small-scale 
magnetic fie l ds has been estab lished from observations 
( KelleJ [l992t [Berger & Title! l200l1). m agneto-hydrody n amica l 

ICarlsson et all [2004; 



simulations 
IShelvag et alj 
conc entrations 
2001a). Other 



(Schiissler et al. 
2004 



2003; 



or semi-emp irical 
dSteiner et all 12001 



models of flux 
ISanchez Almeida et alj 
concentrated on the evolution 



studies have 

of BPs jBerger & Title! Il996t Ivan Ballegooiien etail [l998; 



iBovelet & Wiehi1l2003t iBonet et al. 2005) and on their spectral 
signatures dLanghans et alj|2002ll2004t) . 

Despite these advances, the magnetic field strength and other 
physical properties of the G-band BPs remain largely unknown. 
This is because most studies have used narrow-band filtergrams 
and longitudinal magnetograms, which offer high spatial reso- 
lution but do not allow the determination of the three compo- 
nents of the vector magnetic field. In this paper we investigate 
the magnetic, kinematic, and thermal properties of G-band BPs 
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in the moat of a regular spot using simultaneous vector polari- 
metric measurements of four visible and three infrared lines, 
complemented with diffraction-limited narrow-band filtergrams. 
We also characterize the physical properties of the outer part of 
the moat and surrounding quiet Sun, which do not show G-band 
brightness enhancements. Our aim is to provide a complete ob- 
servational picture of the vector magnetic fields of a large sample 
of BPs to help understand their nature and origin. The results of 
this study may be useful for validating theoretical models and 
MHD simulations of magnetic flux concentrations in the solar 
photosphere, as well as to confirm the mechanisms by which G- 
band bright points are generated. 

The multi-wavelength observations used here have become 
possible only very recently, so we give a detailed account of the 
reduction process and the spatial alignment of the data sets in 
Sect. |2 Our study is based on line parameters extracted from 
the observed Stokes profiles (Sect. 13. 2) , as well as on the in- 
version of the visible and infrared lines (Sect. 13.31 ). In Sect. [4] 
we present the main results of a statistical analysis of the field 
strength, field inclination, magnetic flux, velocity, temperature, 
and size of BPs. Our findings are discussed in Sect. [5] and sum- 
marized in Sect. [6] 

2. Observations and data reduction 

We observed the sunspot NOAA 10425 on August 9, 2003 
from 09:36 to 10:34 UT with the two spectropolarimeters of the 
German Vacuum Tower Telescope (VTT) at the Observatorio 
del T eide (Tenerife, Spaing: the Tenerife Infrared Polarimeter 
(TIP; [Martinez Pifletetaljri 9991) and the Polarimetric Littrow 
Spectrograph (POLIS; ISchmidt et aTJ 120031: iBeck et alj |2005b). 
Both instruments were used simultaneously. The same sunspot 
was observed between 08:25 and 11:58 UT with the Dutch 
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Table 1. Spectral lines observed at the VTT. 



Instrument 


Species 


A a [nm] 




geff 


Pixel size 


TIP 


Fei 


1564.7410 


-0.950 


1.25 


0735 


TIP 


Fei 


1564.8515 


-0.669 


3.00 




TIP 


Fei 


1565.2874 


-0.095 


1.45 




POLIS vis 


Fei 


630.15012 


-0.750 


1.67 


a/ 145 


POLIS vis 


Fei 


630.24936 


-1.236 


2.50 




POLIS vis 


Fei 


630.34600 


-2.550 


1.50 




POLIS vis 


Til 


630.37525 


-1.444 


0.92 




POLIS UV 


Can 


396.849 






07292 



Open Telescope (DOT) at the Observatorio del Roque de los 
Muchachos (La Palma, Spain). The sunspot was located at an 
heliocentric angle of 27°. 

2.1. Observational setup at the VTT 

TIP uses the main spectrograph of the VTT for observations in 
the infrared (IR), while POLIS is a stand-alone spectropolarime- 
ter with two channels. One records the pair of Fei lines at 630 
nm, and the other measures the Ca n H line at 396 nm. In the fol- 
lowing, the two POLIS channels will be referred to as "visible" 
and "UV" channels, respectively. 

An achromatic 50-50 beam splitter was used to split the light 
between TIP and POLIS. Scanning of the solar surface was per- 
formed with the tip-tilt mirror of the correlation tracker sys - 
tem (CT; ISchmidt & Kentischedfl995l iBallesteros etaf] [l996). 
which affects the FOV of both instruments in the same way. The 
CT also provided image stabilization, allowing us to reach a spa- 
tial resolution of about 1", as derived from the power spectrum 
of a granulation area (cf. Appendix [A). Both TIP and POLIS 
measure the polarization state by modulation of the incoming 
light with ferro-electric liquid crystals and a rotating waveplate, 
respectively. To observe the same field of view (FOV) with the 
two instruments, we aligned the TIP and POLIS slits by rotating 
the VTT main spectrograph to the POLIS slit orientation. The 
lateral displacement of the images due to the beam splitter was 
corrected for with the scan mirror of POLIS. 

The slit width was 0735 for TIP and 0748 for POLIS. The 
FOV along the slit was 34" for TIP and 47" for the visible chan- 
nel of POLIS. The step of the spatial scan was 0735, with a total 
integration time of 6 seconds per scan step. The observations 
presented here consist of eight scans of 70 steps (2475) over 
the limb side part of the sunspot and its surroundings. The ca- 
dence for the VTT data is about 7 minutes. Table Q] summarizes 
the properties of the spectral lines observed at the VTT: Ao is 
the laboratory wavelength (taken from Nave et al. 1994 for Fe i) 
and geff the effective Lande factor. The pixel size along the slit 
is given in the last column. The full Stokes vector is available 
for all visible and infrared lines. For Can H, only the intensity 
profile was recorded due to the low light level. FigureQ]displays 
maps of the observed area for the first scan after the coarse align- 
ment described in Appendix |B.2| 

2.2. Observational setup at the DOT 

The data from the DOT consist of speckle-reconstructed filter- 
grams in the G band at 430.5+0.5 nm, Can H at 396.8+0.06 
nm, and blue continuum at 432.0+0.3 nm, with a cadence of 1 
minute. The size of the FOV is 77 "x 60" and the spatial sampling 
07071 per pixel. More detailed information on t he image pro- 
cessing and data reconstruction can be found in ISiitterlin et alj 




Fig. 1. Maps of the first scan across the limb-side part of NOAA 
10425 after coarse alignment of DOT and VTT data. Top, left to 
right: blue continuum, G-band, and Ca n H line-core filtergrams 
from the DOT, and difference between the G-band and blue con- 
tinuum images. Bottom, left to right: continuum intensity at 1.5 
/urn, continuum intensity at 630.4 nm, intensity of the K2V emis- 
sion peak of the Can H line, and intensity in the Can H line 
wing; data from the VTT. 



d2004l) . The angular resolution of the filtergrams reaches the 
diffraction limit of the telescope (072 at 430.5 nm, cf. Appendix 
0). 

2.3. Calibration of polarimetric data 

Besides the usual dark and Hatfield corrections, the polariza- 
tion measurements were corrected for instrumental polarization. 
This was done in two steps. First, TIP and POLIS have inter- 
nal calibration units in front of the CT optics. Using them, the 
crosstalk among the Stokes parameters due to the opt ics behind 
the ca libration units can be determined and removed (Bec k et all 
l2005bl) . The remaining optical path consists of the telescope 
proper, i.e., the two coelostat mirrors, the entrance and exit win- 
dows of the evacuated tube, the main mirror, and one folding 
mirror. The instrumental polarization induc ed by these element s 
was corrected using the telescope model of Bec k et all (|2005a). 
Residual crosstalk between the different Stokes parameters is es- 
timated to be on the order of 10~ 3 / c . The rms noise in the profiles 
is 4 x 10~ 4 / c for the IR spectra and 10~ 3 / c for the visible spectra, 
as evaluated in continuum windows. 

2.4. Alignment of data sets 

For the investigation of the magnetic properties of G-band BPs, 
which are identified using the DOT data, it is crucial to have 
a good alignment of the data sets from the two telescopes. 
Additionally, the simultaneous inversion of infrared and visi- 
ble spectral lines requires the polarization signals in both wave- 
length ranges to be cospatial and cotemporal. Appendix A ex- 
plains in detail the procedures we followed to align the different 
data sets. Figure [2] displays continuum intensity maps from TIP 
and POLIS after coalignment. As can be seen, the spatial cor- 
respondence of the structures in the two images is remarkably 
good. The accuracy of the alignment between TIP and the visi- 
ble channel of POLIS is about 071, or roughly a third of a TIP 
pixel. The maximum displacement perpendicular to the slit due 
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Fig. 2. Continuum intensity maps from POLIS (630 nm; left) and 
TIP (1565 nm; right). Only the lower part of the FOV is shown. 
Contours outline the brightest granules and darkest intergranular 
lanes observed in the POLIS map. Tick marks are separated by 
1". 



to differential refraction (see, e.g.. lReardonll2006h was 1", cor- 
responding to a maximum time gap between cospatial Stokes 
profiles of 18 seconds (3 steps of 0735). 

The time difference between cospatial VTT and DOT obser- 
vations is 30 s or less. The spatial misalignment between VTT 
and DOT data is smaller than one TIP pixel, i.e., OV35. This is 
demonstrated in Fig. [3j where we display coaligned maps for 
the first scan across the FOV. Only the lower part of the FOV 
without the sunspot is shown. The intergranular lanes visible in 
the DOT G-band filtergram can clearly be traced in the IR con- 
tinuum map. Appendix ICl contains all coaligned DOT and VTT 
maps. 

3. Data analysis 

The cospatial and cotemporal observations described in the pre- 
vious section were used to investigate the magnetic, kinematic, 
and thermal properties of G-band BPs in the moat of NOAA 
10425. As a first step, we set up a common wavelength scale 
for the visible and infrared lines. Next, a number of observables 
were extracted from the profiles. These line parameters, together 
with the results of the inversion of the observed spectral lines, 
allowed us to determine the physical properties of the various 
structures in the FOV. Part of our analysis is based on a statistical 
study that compares the properties of BPs and non-BPs (NBPs), 
the latter being points in the FOV that do not show enhanced 
G-band brightness (cf. Sect. 13.41 ). In this section we give details 
on the wavelength scale setup, the extraction of line parameters, 
and the inversion procedure. We also explain how BPs and NBPs 
were identified in the DOT filtergrams and the VTT maps. 

3.1. Wavelength scale 

The wavelength scale depends on two factors: the spectral dis- 
persion and the rest wavelength of the lines, A' Q , where index 
i cycles through the spectral lines observed (cf. Table Q]). The 
dispersion affects, e.g., the determination of the magnetic field 
strength from the separation of the Stokes V peaks, or any ve- 
locity defined by the relative displacement to A' Q . 

To set up a common wavelength scale for TIP and POLIS, we 
calculated a mean quiet Sun (QS) intensity profile in the lower 
half of the FOV. This was done for each of the eight scans and 
each spectral region (IR or visible) separately. We then set the 
rest wavelengths of Fe i 1564.8 nm and Fe 1 630. 15 nm such that 
the observed positions of the line cores in the mean QS pro- 



file were at the predicted convective blueshift values (cf. the last 
row of Table|4]i. Convective blueshifts were calculated using the 
field-f ree, two-component QS model of lBorrero & BellotR ubio 
(2002). Next, the dispersions were adjusted until Fe 1 1565.2 nm 
and Fei 630.25 nm showed the correct convective blueshifts. 
This procedure yielded a dispersion of 29.65 mA pixeL 1 for the 
IR lines and 14.86 mA pixeL 1 for the visible lines. To remove 
the temporal evolution of the relative motions between the Sun 
and the observer, a separate calibration was used for each scan. 

The measurement of Doppler shifts relative to the rest wave- 
length of the lines should be accurate to ±100 m s _1 , as this 
is the uncertainty in determining line-core positions. However, 
the assumed rest wavelengths (and hence the zero points of the 
velocity scales) are more uncertain because of three possible in- 
fluences: a reduction of convective blueshifts due to the presence 
of magnetic fields, a reduction of convective blueshifts because 
of the off-center position of the spot, and a global offset due to 
the moat flow. The first two effects would shift all velocities to- 
wards higher redshifts, because we enforce the lines in the mean 
QS profile to have the convective blueshift values appropriate 
for disc center. However, we believe the error due to a possible 
reduction of convective blueshifts to be small. First, the granula- 
tion appears to be relatively undisturbed in the region where the 
average QS profiles were computed (the moat flow prevents the 
formation of large-scale plage fields close to the spot). Second, 
Balthasar (1985) showed that convective blueshifts outside disk 
center are reduced, but only by a small amount (depending on 
the line). 

With respect to the moat flow, it is important to recall that 
this is a large-scale coherent phenomenon. Since the profiles 
used to compute the average QS profiles are affected by the moat 
flow, referring the velocities to the mean QS profile takes out the 
wavelength shift induced by the moat flow in the same way as all 
other relative motions between observer and target are removed. 
The only influence of the moat flow could then be differential 
moat velocities across the FOV Assuming a ho r izontal moat 
flow of around 0.5 km s _1 (Balt hasar et alJ[T996l : iRezaei et al.l 
2006), we find that the maximum error due to differential ve- 
locities cannot be larger than 0.25 km s at the heliocentric an- 
gle of the observations. The accuracy of the velocity scale thus 
only depends on how appropriate the convective blueshift values 
predicted by the QS model are. As the main effects (off-center 
location, presence of fields) would reduce the blueshift values 
in the QS profiles, all velocities derived from our wavelength 
scale may only be biased to show larger blueshifts than actually 
present. 

We note that, prior to the determination of the wavelength 
scale, the profiles from the visible channel of POLIS have been 
corrected for the curvature of the spectrum along the slit. The 
curvature is caused by the short focal length of the spectrograph 
and the large FOV along the slit. To avoid errors in the velocity 
determination due to this effect, the spectra have been shifted in- 
dividually so as to have the teluric O2 630.20 nm line (which 
does not show Doppler velocities) at the same position as in 
the QS profile from w hich the wavelength scale is determined 
(cf. IRezaei et alj|2006l) . 

3.2. Line parameters 

For each spatial point, we have the four Stokes profiles of the 
visible and IR lines listed in Table [1] the intensity profile of the 
Ca n H line from the POLIS U V channel, and the intensity in the 
G Band, the Ca n H line core, and the blue continuum from the 
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Fig. 3. Left panels: G-band intensity (top) and IR continuum intensity (bottom) for the first scan. The intergranular lanes visible in 
the G-band map can be traced in the IR intensity map as well (black contours) despite the inferior spatial resolution. Tick marks are 
separated by 1". Right panels: Total integrated polarization (top) and DOT BP mask (bottom). The individual BPs are color coded 
in the DOT BP mask. White contours outline strong polarization signals (> 1%) in the moat of the spot, as well as the sunspot 
canopy. Most of the identified BPs (red contours) coincide with areas of non-zero polarization signal (top right). 



DOT, averaged over a 5 x 5 pixel area to simulate a TIP pixel. 
After alignment of the data sets, these observables are cospatial 
and cotemporal. 

The four spectral lines with the highest magnetic sensitivity 
(Fei 1564.8 nm and Fei 1565.2 nm in the IR, Fei 630.15 nm 
and Fe i 630.25 nm in the visible) have been used to extract the 
following line parameters: 

- The line-core velocity from the intensity profiles (also for 
Tii 630.37 nm). 

- The total integrated polarization, T, defined as 



T 



(Q £ + U 2 + V 2 ) l/2 /I c dA 



(1) 



in a wavelength range encompassing the full spectral line. 
- The maximum polarization degree, 



p = max{(e 2 + t/ 2 + V 2 ) 1/2 //}- 



(2) 



A profile is labeled unpolarized if the maximum polariza- 
tion degree is below 0.2% (IR lines) or 0.5% (visible lines), 
otherwise it is called polarized. These thresholds correspond 



to 5 times the noise level. Pixels with at least one spectral 
line above the respective threshold were subject to a two- 
component inversion (cf. Sect. |3.3| below). 

The following line parameters were calculated only for profiles 
showing polarization signal: 

- The mean ratio of total linear and circular polarization, 
L/V = ((Q 2 + U 2 ) 1 ^ 2 /\V\), in a narrow wavelength range 
around the maximum absolute Stokes V signal, i.e., the ratio 
at the position of the strongest <x component. 

- The area asymmetry of the Stokes V profile, 



6 A = s x 



J V(A)dA 
}\V(A)\dA 



(3) 



The polarity, s, is ±1 and denotes the orientation of the mag- 
netic field vector to the LOS. It is derived from the position 
of minimum and maximum Stokes V signal. 
The number and position (with pixel accuracy) of lobes in 
Stokes Q, U, and V. Lobes are defined as local extrema of 
smoothed polarization profiles above a threshold of 0. 1 % for 
the IR lines and 0.25% for the visible lines. 
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Fig. 4. Example of cospatial Stokes profiles from a G-band BP, normalized to the continuum intensity of the quiet Sun at disk center. 
From top to bottom: Stokes /, Q, U, and V. The vertical lines mark the rest wavelengths derived from the quiet Sun intensity profiles. 
The observed profiles are indicated by filled circles; the solid lines represent the best-fit profiles from the inversion. Left panel: Fe i 
line at 1564.8 nm with blend at 1564.7 nm, and Fe i 1565.2 nm. Right panel: Fe i lines at 630. 15 and 630.25 nm with the teluric O2 
lines at 630.200 and 630.276 nm, Fe 1 630.35 nm and Ti 1 630.37 nm. The Stokes / and V profiles of the Ti 1 line are magnified by a 
factor of 5 for better visibility. 



Finally, for each Stokes V profile with two distinct lobes, we 
calculated 

- The Stokes V amplitude asymmetry, 

Iflbluel - l^redl 

5a = — r, (4) 

Iflbluel + lOredl 

where flbiue and a re( j represent the amplitudes of the blue and 
red lobes, respectively. 

- The wavelength position of the V lobes with sub-pixel accu- 
racy for the derivation of field strengths. 

- The Stokes V zero-crossing velocity, v zcro , from a linear fit 
to the circular polarization signal in between the lobes. 

3.3. Inversion procedure 

The four visible and the three infrared lines observed at the VTT 
were inverted together using the SIR code (Stokes Inversion 
based on Response functions; Ruiz Cob o & del Toro Iniestal 
119921) . To facilitate the recovery of temperatures in deep photo- 
spheric layers, the Stokes profiles of each line were normalized 
to the continuum intensity of the quiet Sun at disk center. In ad- 
dition, we removed the teluric blend of the Fe 1 1564.8 nm line in 
the following way. The Fe 1 1564.8 nm quiet Sun intensity profile 



was inverted with a two-component quiet Sun model. The best- 
fit profile does not contain the teluric blend, but reproduces the 
continuum outside the spectral line. The division of the observed 
with the best-fit profiles yields a multiplicative correction for the 
removal of the blend that was applied to all Fe 1 1564.8 nm pro- 
files. 

The model atmosphere and corresponding free parameters 
depended on the type of pixel inverted (Tabled. The free param- 
eters in the case of field-free atmospheres were the temperature, 
T, and line-of-sight (LOS) velocity, v. For magnetic atmospheres 
we also had the field strength, B, the LOS inclination, y, and the 
azimuth of the field in the plane perpendicular to the LOS, i/f. The 
model atmosphere was given as a function of continuum optical 
depth, r, in the range from log r = 1 to -4. The temperature was 
determined using two nodes, i.e., the absolute value and slope of 
the initial stratification could be modified. All other atmospheric 
parameters were assumed to be constant with depth. Our choice 
of height-independent magnetic and kinematic parameters has 
been made on purpose to facilitate the analysis of such a large 
data set. Model atmospheres with height-independent parame- 
ters cannot reproduce the asymmetries of the observed Stokes 
profiles, but have been shown to provide reasonable averages of 
atmospheric parameters along the LOS. This is exactly what is 
needed for a statistical analysis like the present one, where we 
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Table 2. Free parameters for the inversion of unpolarized pixels 
(second column) and polarized pixels (third and fourth columns). 
The meaning of the symbols is as follows: x, free parameter; -, 
not used; 0, same value in both atmospheric components. 



Parameter 


Unpolarized 


Polarized 


Polarized 






1st comp 


2nd comp 


Temperature (T) 


X 


X 


X 


Velocity (v) 


X 


X 


X 


Field strength (B) 






X 


Field inclination (y) 






X 


Field azimuth (if/) 






X 


Microturbulence (Vmi C ) 


X 


X 


X 


Macroturbulence (v mac ) 


X 








Straylight factor (J3) 










Filling factor (/) 











Fig. 5. Atmospheric parameters deduced from the inversion of 
the seventh map. Top: G-band intensity and IR continuum inten- 
sity. Middle: LOS velocity of the magnetic component and mag- 
netic filling factor. Bottom: Field strength and field inclination to 
the LOS. White contours outline areas of enhanced polarization 
signal and red contours (green in the velocity map) the identified 
BPs. Tickmarks are arcsec. Black areas represent pixels with- 
out polarization signal that have not been inverted in terms of a 
two-component model. 

focus on the general properties of G-band BPs rather than on the 
details of line formation. 

The synthetic profiles were convolved with a macroturbulent 
velocity v mac . In our case, v mac describes the effects of large- 
scale motions and the unknown point spread functions (PSFs) 
of the spectrographs. This treatment, however, is simplistic be- 
cause TIP and POLIS have slightly different PSFs. The micro- 
turbulence, v m i c , was another free parameter of the inversion. 
The code determines v mac and v m ; c separately for each inverted 
pixel. However, when two-component model atmospheres were 
used, the value of v mac was forced to be the same in both com- 
ponents. We employed one-component, field-free model atmo- 
spheres to invert pixels with no polarized profiles. The initial 
temperature stratification was t hat of the Harvard Sm ithsonian 
Reference Atmosphere (HSRA. lGingerich et alJH971[) . No stray 
light was allowed for. The free parameters of the inversion were 
T, v, v mac , and v mic . 

Pixels that were classified as polarized (cf. Sect. 13.2b were 
inverted with a two-component model. The first component was 
identical to the field-free atmosphere described above. The sec- 
ond component was magnetic. The Stokes profiles emerging 
from the two components, p m - dg and p nmag , are combined to yield 

PmodelW = (1 - f) Pnmag(T) + / PmagW, (5) 

where / represents the magnetic filling factor, i.e., the frac- 
tional area of the pixel occupied by the magnetic component. 
Additionally, some amount of stray light contamination was as- 
sumed. We adopted the quiet Sun intensity profile as an approxi- 
mation to the stray light profile. In this way, the synthetic profile 
used to fit the observations is 

p(X) = B p Amy (X) + (1 - 8) Pmode i(A), (6) 

where B represents the stray light factor. This combination is 
intended to describe Stokes profiles emerging from unresolved 
flux concentrations and their field-free surroundings, affected by 
stray light coming from unresolved granules and intergranular 



lanes. It is important to point out that B and / have different ef- 
fects on the position and shape of the intensity profiles, which 
makes it possible to distinguish between them. The stray light 
profile corresponds to an additional field-free atmospheric com- 
ponent with fixed velocity and fixed temperature stratification. 
In contrast, the field-free component of the inversion can be ad- 
justed both in velocity and temperature to represent the physical 
conditions of the plasma in the immediate vicinity of the flux 
concentration. We have chosen the simplest model possible for 
the inversion in order to have a robust estimate of magnetic prop- 
erties of G-band BPs, which of course has some drawbacks. For 
example, the restriction to parameters that are constant in opti- 
cal depth does not allow to reproduce asymmetries in the Stokes 
profiles. 

The uncertainties in the atmospheric parameters resulting 
from the inversion were es timated through th e diagonal elements 
of the covariance matrix dPress et al.| [T986). Averaged over the 
4000 inverted profiles of the first scan, the formal errors in B, 
y, v, and / turned out to be 100 G, 10°, 400 m s" 1 , and 2%. 
A number of factors may contribute to the relatively large error 
in velocity: possible height-variations of the flows that cannot 
be accounted for by our simple model, errors in the convective 
blueshift values (+150 m s _I for each line), and uncertainties in 
the laboratory wavelengths (especially for Tii 630.37 nm). Yet, 
these errors are smaller than the uncertainties associated with 
magnetogram or Dopplergram observations. The reason is the 
large amount of information we have at our disposal: the four 
Stokes profiles of seven different lines, which are used to con- 
strain the model parameters simultaneously. Examples of ob- 
served and best-fit profiles are displayed in Fig. [4] Maps of the 
atmospheric parameters resulting from the inversion are shown 
in Fig.|5] 

3.4. Identification of BPs 

In order to identify BPs in the DOT data, we calculated the rela- 
tive G-band enhancement or G-band contrast map, C, using the 
G-band intensities, /gb, and blue continuum intensities, /bc, as 
follows: 

r _ Igb • (he) /{Igb) - he 
Vbc) 

where (/) denotes the spatial average of the intensity over the 
map. Only strong enhancements of the G-band intensity relative 
to the blue continuum intensity are left over. An example can be 
found in Fig.Q] 

The BPs proper were taken to be all closed contours in the 
G-band contrast map above a threshold of 0.31. This value was 
selected by trial and error so as to encompass all BPs an observer 
would have identified visually in the G-band images of the first 
scan. Our procedure misses G-band BPs with low contrasts, but 
the detection of these BPs is always problematic independent of 
the method used to identify them. Contours near the upper edge 
of the FOV were not considered in the analysis to avoid contam- 
ination by the sunspot canopy. A few very short contour lines 
extending over 1-5 DOT pixels were excluded manually since 
they are presumably due to spurious signals created by the align- 
ment procedure. Contours tracing G-band brightenings inside or 
at the edge of granules were removed as well. No further crite- 
rion on the contrast or the shape of the areas was applied. The 
closed contours were consecutively numbered to address indi- 
vidual BPs. In total, 447 separate BPs were identified in the full 
data set. Our selection of bright areas by sharp thresholding leads 
to extended patches, which would perhaps be split into smaller 
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Fig. 6. Top: Relative frequency of total integrated polarization 
values for BPs (thick black) and NBPs (thin grey). Middle and 
bottom: G-band contrast vs total integrated polarization for BPs 
and NBPs. The dashed horizontal lines indicate the mean con- 
trast of 0.52 for BPs and 0.01 for NBPs. The shaded areas denote 
the range of T values within the noise level. 



patches by t he use of more sophis ticated methods like the one 
proposed bv lBovelet & WiehJ J2003I) . 

A BP mask for the polarimetric data from the VTT was con- 
structed by marking all VTT pixels cospatial to one or more 
pixels above the contrast threshold in the DOT data. The 447 
BPs identified in the G-band contrast maps extend over 1238 
VTT pixels. The control sample of non-bright points (NBPs) was 
taken to be all VTT pixels at least three pixels (l'/05) apart from 
any BP and outside the canopy of the spot (visible near the upper 
edge of the FOV). An example of the NPB mask is displayed in 
Fig. lC.ll for the first map. In total, around 9500 pixels were se- 
lected for the NBP sample. 

The majority of BPs identified by the DOT mask have sizes 
below the spatial resolution of the polarimetric data. Part of 
the polarized light produced by these unresolved sources will 
therefore be distributed over several VTT pixels according to 
the spatial point spread function. The construction of the mask 
for the VTT data thus corresponds to choosing the spectra of 
the location where the signature of the BP should be maximum. 
However, light from other polarization sources inside the reso- 
lution element not related to a specific BP can also contaminate 
the spectra attributed to the BP. 



Table 3. Fraction of unpolarized and polarized profiles of differ- 
ent types associated with BPs and NBPs, in percent. 

BP profiles, 1238 in total 



Line 


Unpolarized 


Polarized 




(nm) 




Total Regular Irregular Blue only Red only 


1564.8 


3 


97 (35) (61) < 1 


< 1 


1565.2 


11 


89 70 16 3 


<1 


630.15 


12 


88 80 7 < 1 


1 


630.25 


5 


95 83 6 3 


3 


NBP profiles, 9497 in total 


Line 


Unpolarized 


Polarized 




(nm) 




Total Regular Irregular Blue only Red only 


1564.8 


33 


67 (33) (26) 6 


2 


1565.2 


72 


28 15 4 7 


2 


630.15 


74 


26 20 < 1 3 


2 


630.25 


49 


51 22 2 20 


7 



3.5. Data representation 

In the next sections, the parameters we use to describe the indi- 
vidual BPs are averaged over the spatial extension of each BP in 
the DOT mask and, analogously, over the corresponding area in 
the VTT mask if the BP extends over more than one pixel. The 
NBPs always encompass single pixels, so no averaging is carried 
out. 

Individual G-band BPs may show different physical proper- 
ties (e.g., magnetic flux, field strength, LOS velocity); in addi- 
tion, there is noise in the data. Both factors lead to an unavoid- 
able scatter that hides possible relations between the parameters 
characterizing the BPs. In an attempt to reveal these relations 
more clearly, we have chosen to investigate mean parameters 
obtained by binning the values of individual BPs. For all scat- 
ter plots relating two quantities, e.g., field strength and G-band 
contrast, we sorted the pairs (x,y) of abscissa, x, and ordinate, y, 
to be in ascending order in x. We then averaged a variable num- 
ber, N, of values in x (typically N — 5 - 50), and averaged the 
values of the dependent variable y over the same pairs. With this 
procedure we keep the same statistics for each point plotted, as 
it is always derived from the same number of (x,y)-pairs. The 
range in x over which the averaging is performed depends on 
how dense the values in x are distributed after sorting. 

Thus, each data point represents an average, and the error 
bar indicates the uncertainty of the mean, crj yN, where cr is the 
standard deviation of the averaged values in y. The implicit as- 
sumption here is that for a trend of, e.g., G-band contrast with 
field strength, other quantities like field inclination are randomly 
distributed, and their influence cancels out by the averaging pro- 
cess such that only the dependence on field strength is left over. 
Whenever possible, the noise level is estimated from the sample 
of NBPs. In most cases, the corresponding results for the NPBs 
are also shown for comparison. The statistics of the NPBs are 
better due to the larger number of data points. 



4. Results 

Our results derive either directly from the spectral line profiles 
or from the more complex inversion procedure. In some cases, 
we used both approaches to confirm specific findings, while in 
other cases they were complementary. 
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Fig. 7. Magnetic field strength of BPs from the Zeeman splitting 
of Fei 1564.8 nm (+), 1565.2 nm (□), and 630.25 nm (A), vs 
the values inferred from the inversion. The shaded areas give the 
scatter inside the bins. The straight line indicates a one-to-one 
correspondence. 



4.1. Polarization signature of BPs 

A visual inspection of Figs.[3][5] or lC.ll suggests that the majority 
of G-band BPs are associated with polarization signals. To quan- 
tify this impression, we use the total integrated polarization, T . 
T is calculated separately for the four most Zeeman-sensitive 
spectral lines according to Eq. ([TJ, and then averaged over all 
lines. We estimate the noise level in T by averaging T over all 
profiles that were not inverted due to their maximum polarization 
degrees being below the corresponding threshold. The average 
value plus one standard deviation, T mem +<T, is 0.22 pm for these 
profiles. Values below this limit are assumed to correspond to the 
integration of pure noise, where of course we cannot exclude the 
existence of weak polarization signals below the noise level. The 
histogram of T values (upper panel of Fig. |6]l demonstrates that 
the majority of BPs (94%) have polarization signals above this 
level, whereas half of the NBPs have T values below it. Thus, 
clear polarization signals (i.e., magnetic fields) are observed no 
farther than ff.'35 from the position of almost all G-band bright- 
enings. Due to the 1 "spatial resolution of the spectra we cannot 
exclude that the polarization signal comes from a source outside 
the pixel itself, but we think it highly improbable that it happens 
so often by chance. The middle and lower panels of Fig. [^dis- 
play the G-band contrast of BPs and NBPs as a function of T. 
We find that the contrast of BPs only depends on T for small 
polarization signals on the order of the integrated noise. For T 
values higher than 0.3 pm, the BP G-band contrast is approxi- 
mately constant. The G-band contrast stays close to zero for the 
NBPs and shows no systematic trend with T . This suggests that 
the NBP sample is distinctly different from the BP sample, even 
if it contains locations that exhibit some polarization signal. 

The shape of the Stokes V profiles gives information about 
the geometry of the magnetic field and the dynamical state of the 
flux concentrations. After the distinction between polarized and 
unpolarized profiles, the polarized profiles (cf. Sect. 13.21 ) were 
classified according to the number of Stokes V lobes. A profile 
is called regular if it has 2 lobes, irregular if it shows 3 or more 
lobes, and blue/red only if it shows only one blue/red lobe at 
shorter/longer wavelength than the intensity minimum. The Fe i 
line at 1564.8 nm is blended with Fe 1 1564.7 nm, which explains 
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Fig. 8. Top: Distribution of BP field strengths as inferred from 
the inversion (thick line) and the splitting of the Fei 1564.8 
nm line (thin line). Middle: same for the NBP sample. Bottom: 
Fraction of the total area occupied by fields of given strength. 



the high number of irregular profiles; the respective numbers are 
less reliable and thus are given in parentheses (Table|3]l. We find 
that roughly 80% of the BP profiles are regular with two lobes in 
Stokes V, while 10 % are irregular with either one or three lobes, 
or unpolarized (10%). In the sample of NBPs, the fraction of un- 
polarized and irregular Stokes V profiles is significantly larger. 
A high percentage of BPs are associated with regular two-lobed 
Stokes V profiles. There are few cases of irregular three-lobed 
profiles, which would indicate fields with opposite polarity and 
different LOS velocities in the same pixel. Examples of this kind 
can be found in the sunspot penumbra at the upper boundary of 
the FOV. Thus, even if the asymmetries of the Stokes V profiles 
(cf. Sect. 14. 5t indicate some variation in the magnetic properties 
along the LOS, most BPs do not seem to be related to a complex 
multi-component field topology. 



4.2. Magnetic field strength 

For spectral lines in the strong field limit, where the separation 
of the cr components is larger than the thermal Doppler broad- 
ening, the magnetic field strength (B) can be calculated from the 
position of the Stokes V lobes according to 



AA = 4.67 x lQ- u g eS BAl 



(8) 



where AA is the wavelength distance between the peaks of Stokes 
V and the line center. In observations of magnetic fields outside 
sunspots or pores, Eq. ([8]) only fully applies to the IR line at 
1564.8 nm for fields above ~500 G. However, for a comparison 
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Fig. 9. G-band contrast vs field strength for the BP sample. The 
dashed line represents a linear fit to the data points. 



with the inversion results, we calculate B from Eq. ([8]l for the 
three most Zeeman-sensitive lines. 

Figure|7]displays the BP field strengths indicated by the split- 
ting of the different lines versus the field strengths inferred from 
the inversion. Taking the latter as the most accurate ones, it is 
clear that the splitting of Fei 1564.8 nm (+) serves as a good 
diagnostics of B for fields down to 500 G. The Fei 1565.2 nm 
line (□), and the visible line at 630.25 nm (A) especially, show 
strong deviations for fields below 1200 G. The same result was 
obtained for the NBP sample. 

The distribution of field strengths is displayed in Fig. [8] for 
individual BPs and NBPs. The inversion and the splitting of Fe i 
1564.8 nm essentially give the same result: the field strengths 
of BPs are distributed uniformly in the range from 500 G to 
about 1300-1400 G. The finding of G-band BPs with strengths 
of ~ 500 G may pose a serious problem for our understanding 
of why these flux concentrations are bright. In general, mod- 
els predict zero or even negat i ve G-band contrast for 500 G 
dSanchez Almeida et al. 2001a; Shelv ag et all |2004|) . although 
high positive contrasts are also detected. The distribution of field 
strengths in the NBP control sample hints at a quickly decreasing 
proba bility for stronger fie lds. The same result has been found 
by iKhomenko et aTl d2003l) in their study of quiet Sun internet- 
work fields using IR lines. The sharp drop of the relative fre- 
quency of fields weaker than 300 G may be imposed by th e 
detection limit of the observations (cf. IKhomenko et alj |2003). 
Remarkably, it seems that the spot does not influence the distri- 
bution of weak fields in the outer part of the moat and surround- 
ing quiet Sun much (the regions where the NBP sample comes 
from, cf. Fig. IC. 11 1. To allow comparisons with magnetoconvec- 
tion simulations, the lower panel of Fig.|8]shows the area occu- 
pied by magnetic fields of a given strength using 100 G bins, as 
a fraction of the total area occupied by fields in the BP and NBP 
sample, respectively. Figure [9] demonstrates that there is a slight 
increase in the G-band contrast toward higher field strengths for 
the BPs. 

4.3. Magnetic field inclinations 

Figure [10] displays histograms of the LOS inclination for BPs 
and NBPs. Positive and negative polarities are found in both 
types of structures. In the case of the NBPs, the inclinations 
are distributed rather symmetrically around two broad peaks at 
y = 40° and 130°, with few fields perpendicular to the line of 
sight. For comparison, we also show the histogram of LOS in- 
clinations found at the upper border of the FOV, close to the 
penumbra and its filaments (cf. Fig. IC.lt . The fields in this area 
exhibit a single peak at around 70°, which corresponds to the 
unipolar horizontal fields of the sunspot canopy. The distribu- 



tion of LOS inclinations for the BPs peaks at around 15° and 
145°; the peaks are less broad than those of the NBP sample. 
Few BPs show fields perpendicular to the LOS. The LOS incli- 
nations of BPs and NBPs at a heliocentric angle of 9 — 27° are 
consistent to the first order with fields that are close to vertical 
to the surface. The systematic displacements toward smaller in- 
clination (BP peak at 15°, NBP at 130°) could be due either to 
the influence of the c anopy fields of the sunspot or to unresolved 
mixed polarity fields (Lites 2 0021) . 

The inclination of the vector magnetic field to the LOS can 
also be estimated from the ratio of linear to circular polarization, 
L/V. The L/V ratio is more affected by noise in the data than 
the inclination derived from the inversio n, especially for small 
amplit udes of L or V. For fully split lines, lLandi DeglTnnocentil 
(2003) demonstrate that 



1 sin y 

2 cosy 



(9) 



The upper left panel of Fig. QT]displays a scatter plot of the 
L/V ratio for the two most Zeeman-sensitive spectral lines, Fe i 
1564.8 nm and Fe i 630.25 nm. The values of L/V are similar in 
both wavelength regions; they only show significant deviations 
for values of L/V above 0.45 or, according to Eq. (0, for inclina- 
tions larger than 50°. The upper right panel of Fig. [TTIcompares 
the L/V values with the results of the inversion. This plot demon- 
strates that Eq. (0 is also valid for the visible Fe 1 630.25 nm line 
if the field is not too much inclined (y <k 50°). In general, the 
observed values of L/V are slightly higher than those indicated 
by the inversion, which we attribute to the influence of noise. 

Interestingly, the G-band contrast of BPs seems to depend on 
the magnetic field inclination. The middle panels of Fig.fTTIshow 
the contrast as a function of both L/V and the magnetic field in- 
clination inferred from the inversion. We find a correlation of 
the brightest structures with magnetic fields parallel to the LOS. 
The BPs with magnetic fields perpendicular to the LOS system- 
atically exhibit reduced G-band contrasts. We believe this trend 
to be significant even if the contrast shows some scatter due to in- 
trinsic differences in the properties of the BPs like field strength 
or flux, which also have some effect on the contrast. The over- 
plotted curves will be discussed in Sect. [5] No relation between 




Fig. 10. Histogram of the magnetic field inclination with respect 
to the LOS for BPs (black), NBPs (grey), and the canopy area 
close to the sunspot (light grey). The vertical dashed line cor- 
responds to fields perpendicular to the LOS. The vertical solid 
lines mark the LOS inclinations of fields, which are purely ver- 
tical relative to the solar surface. 
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for three values of the parameter b. Error bars indicate the uncertainty of the mean values. Bottom panels: BP contrast vs L/V and 
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the G-band contrast and the field inclination is observed for the 
NBPs (bottom panels of Fig.Qj}- 

4.4. LOS velocities 

We used two proxies to estimate flows of magnetized and field- 
free plasma: the Stokes V zero-crossing shift and the Stokes / 
line-core velocity, respectively. The inversion also yielded values 
for the LOS velocity in the magnetic and field-free components 
of the atmosphere. By convention, positive velocities indicate 
redshifts. 

Figure [12] shows the dependence of the LOS velocity on the 
G-band intensity. For easier comparison with other data sets, we 
used the relative G-band intensity 1 + C in this plot. For the 
NBP sample (small symbols in the left half of the figure, with 
relative intensities below 1.31), we find that the velocities de- 
rived from the Stokes / line-core position and the velocities of 
the field-free component of the inversion are in good agreement 
with the expectation for an area dominated by granulation: bright 
granules show blueshifts (v < 0), whereas intergranular lanes 
have reduced intensity and are associated with redshifts (v > 0). 



In contrast, no trend of the magnetic velocity with intensity is 
found in the BPs (large symbols, relative intensities above 1.31 
in Fig.[T2l. The velocity of the BPs spans the range ±2 km s" 1 for 
the binned data. The extreme velocities, considering only points 
with a clear polarization signal above 1%, reach ±4 km s _1 . In 
general, all spectral lines and the inversion show a preference for 
downflows of magnetized plasma inside the BPs. 

The magnetic velocity indicated by the inversion usually 
agrees with the Stokes V zero-crossing velocity deduced from 
the visible lines (upper panel of Fig.[T3l. because a displacement 
between the observed and synthetic Stokes V profiles would 
strongly degrade the quality of the fit, due to the steep slope of 
the circular polarization signal of the visible lines near the zero 



Table 4. LOS velocities of BPs and NPBs in m s" 1 . 
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1564.8 


1565.2 


BPs, zero-crossing 
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NBPs, line core 
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Convective blueshift 


-185 


-262 


-424 


-445 


-469 



C. Beck et al,: Magnetic properties of G-band bright points 



11 



0.4 
0.2 
0.0 
-0.2 
-0.4 











c^ + * ++ ++ + + ++ + + 


+ 

4- 


+ + 

3><X> 

'■ + 


+ 






o 1565.2 " 
* 1564.8 - 
+ 630.25 ^ 







0.4 
0.2 
0.0 
-0.2 
-0.4 



V amplitude [%] 



2 3 
V amplitude [%] 



o * a 




* o * 


o : 



0.4 

0.2 

0.0 
-0.2 
-0.4 

-3000 -2000 -1000 



1000 2000 3000 
, [m/s] 



o 


o 


: o + » + +" e 





0.4 

0.2 

0.0 
-0.2 
-0.4 

-3000 -2000 -1000 



1000 2000 3000 
, [m/s] 




Ampl. asymmetry V amplitude [%] 

Fig. 15. Scatter plot of various Stokes V line parameters for the BP sample. Left to right, top to bottom: 6A and da vs Stokes V 
amplitude, SA and 5a vs zero-crossing velocity, area asymmetry vs amplitude asymmetry, and zero-crossing velocity vs Stokes V 
amplitude. 





2000 
1000 


- t *' 6'36.i5 
-□□ 630.25 
i 630.37 
1564.8 
;.. 1565.2 
7 Inversion 


' ' ' ' ' ' ' □ ' ' 11 ' ' 

□ 

+ A Q 

o f? n + 


□ 







: I 


+ ; 


velocity 




+ □ ,0 

h + n 

fe + + iP ° 
^ Q □ 






-1000 




□ 






-2000 




, , , 1 , , , 1 , ., , 1 , 4. , 1 





0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
rel. G — Band intensity 

Fig. 12. Relation between the LOS velocity and the normalized 
G-band intensity. Positive velocities indicate redshifts. Small 
symbols: Stokes / line-core velocity and LOS velocity of the 
field-free component of the inversion for the NBP sample. Large 
symbols: Stokes V zero-crossing velocity and LOS velocity of 
the magnetic component of the inversion for BPs. 



crossing. The IR line at 1565.2 nm usually shows the strongest 
redshifts. Fei 1564.8 nm has less reliable statistics, as v zclo is 
only calculated for regular profiles (Table [3]). The bottom panel 
of Fig. [13] compares the Stokes V zero-crossing velocities de- 
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Fig. 13. Top: Stokes V zero-crossing velocity, v zcro , vs the mag- 
netic velocity from the inversion, for BPs. Bottom: BP zero- 
crossing velocities deduced from IR lines vs those from visible 
lines. The zero-crossing velocity is averaged over the two visi- 
ble (or infrared) lines when both show regular Stokes V profiles. 
The shaded area indicates the scatter in the bins. 



rived from the IR and visible lines. The plot also shows a trend 
toward higher velocities in the IR as compared with the visible, 
as most points lie above the line of one-to-one correlation. This 
difference in magnetic velocity could be due to unresolved struc- 
tures or flow fields, whose effect on infrared and visible lines 
may differ. It is also compatible with the existence of velocity 
gradients along the line of sight, as indicated by the profile asym- 
metries discussed in the next section. 
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Fig. 14. Histogram of the LOS velocity of the magnetic compo- 
nent of the inversion, for BPs. 

Table|4]summarizes the average velocities of BPs and NBPs; 
positive values correspond to redshifts. The first and second row 
show the average zero-crossing velocity and the average line- 
core velocity of BPs for the different spectral lines. The third row 
contains the average line-core velocity of NPBs. In the fourth 
row, the con vective blueshifts from the tw o-component quiet 
Sun model of lBorrero & Bellot Rubiol (2002) are given that were 
used to set up the wavelength scale. For BPs, the average zero- 
crossing velocities are between 300 m s and 600 m s , de- 
pending on the spectral line; the average velocity of the magnetic 
component of the inversion is 260 m s" 1 (cf. Fig.[T4jl. Th ese ve- 
locities agree with the results of lAmer & Kneeri (1 19931) . as far 
as different velocities are deduced from St oke s / and Stokes 
V, and with iGrossmann-Doerth et al.l d 19961) or ISigwarfh et al] 
(1999), who found significant redshifts in magnetic elements. 
The line-core velocity of the BPs shows blueshifts. The mag- 
netic elements and intergranular lanes are not resolved in the 
VTT data; thus, the intensity profiles of the BPs are probably 
affected by unpolarized light from bright granules in the im- 
mediate surroundings, which is a source of blueshifts. For the 
NBP sample, we find average velocities close to the convective 
blueshifts used in the determination of the wavelength scale, as 
expected. As discussed in Sect. 13.11 the velocities may be biased 
towards the blue. 



and internetwork areas. Negative asymmetries are encountered 
only for weak polarization signals. Both 6A and 8a change sign 
with the zero-crossing velocity (middle panels of Fig. [TBI . As 
a consequence, there are very few examples of, e.g., positive 
asymmetries associated with blueshifted profiles. The lower left 
panel of Fig.[l5]shows that the correlation between area and am- 
plitude asymmetry is quite tight for G-band BPs: structures with 
positive area asymmetry also have positive amplitude asymme- 
try, and when one increases the other follows it closely. Another 
finding is that the zero-crossing velocity of BPs depends on the 
amplitude of the polarization signal only for very low signals: 
most BPs with strong Stokes V signals show similar velocities 
(cf. the lower right panel of Fig. [T3J. 

4.6. Magnetic flux 

The magnetic flux per pixel is calculated from the inversion re- 
sults as O = (1 - P)f A B cosy, where A - 254 2 km 2 is the 
area corresponding to a VTT pixel. Figure [T6]displays the vari- 
ation of the G-band contrast of BPs with the average unsigned 
flux per pixel, (€>)bp, where ()bp indicates the average over all 
VTT pixels covered by the BP. The contrast increases with mag- 
netic flux until about (<1>)bp = 5 x 10 16 Mx and then decreases. 
The relation between G-band contrast and mean unsigned flux 
is very similar to that between contrast and total integrated po- 
larization (T, cf. Fig. [6). This supports the consistency of the 
inversion results: T is proportional to the flux and is derived di- 
rectly from the observed profiles. The magnetic flux is calculated 
using three different free parameters of the inversion, but shows 
the same dependence with the G-band contrast. 

The total magnetic fluxes we infer for the BPs, <J> tot = 
n ■ (O)bp, with n the number of V TT pixels cover ed by the BP, 
are similar to the ones obtained bv lBerger & TMel d200ll) . albeit 
generally smaller (see left panel of Fig. IT7T ). For a comparison 
with their paper we also computed the maximum value of the 
flux density in each BP (O max (BP)/A; right panel of Fig. [TTT i. for 
which we again find lower values than they do. The most prob- 



4.5. Area and amplitude asymmetries 

The area and amplitude asymmetries are sensitiv e to velocity and 
magn etic field gradients along the line of sight. lAuer & Heaslevl 
d 19781) demonstrated that the existence of velocity gradients 
along the LOS is a necessary and sufficie nt condition for hav - 
ing non-zero area asymmetries (cf. also Lope z Aristd [2002). 
Enhanced asymmetries result when gradients of velocity are 
combined with gradients of magnetic field strength, field inclina- 
tion, and/or field azimuth. Note that jumps of these atmospheric 
parameters at the interface between magnetic flux concentra- 
tions and their field-free surroundings would effectively produce 
gradients along th e LOS . I n the context of magnet i c cano pies, 
Solanki & Pahlke (1988}), ISanchez Almeida et al.1 dl989l) . or 
Sanchez Almeida ( 1998) have demonstrated that the sign of the 
area asymmetry is related to the sign of the gradients of velocity 
and field strength through 



sgn (6A) 



/dB(r)dv(T) 

sgn 

6 dr dr 



(10) 



In Fig. [15] we display the correlation between the asymme- 
tries and other parameters for the BPs. The upper panels show 
that most BPs exhibit positive area and am plitude asymmetries, 
with v alues comparable to th ose found by | Martmez Pillet et al.l 
d 1997b in plage regions and by Sig warth et al.l d 19991) in network 




Fig. 16. G-band contrast vs average single-pixel magnetic flux 
of BPs. The curve has a similar shape to the contrast vs total 
integrated polarization (cf. Fig.[6]i. 
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Fig. 17. Left: Histogram of the total unsigned flux of BPs. Right: 
Maximum flux density of BPs. 
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Fig. 18. Top left: G-band contrast vs temperature of the two atmospheric components used to invert the BP profiles: the magnetic 
interior (O) and the external, field-free surroundings (+). Top right: Histogram of temperatures. The magnetic interior is usually 
hotter than the surroundings by about 1000 K at the same optical depth. Bottom left: Temperature of the two atmospheric components 
vs average unsigned magnetic flux. The magnetic temperature is seen to decrease linearly with flux. Bottom right: Temperature as 
a function of field strength. A slight decrease in the magnetic temperature is observed. For the field-free component, no clear trends 
with <5 mean or B exist. 



able value of (D max (BP)/A is around +60 G for positive-polarity 
BPs and -30 G fo r negative-polar i ty BPs , i.e., only about half the 
value reported by Berge r~& Title! d200ll) . This systematic differ- 
ence may have two origins: (1) the different methods used to de- 
rive the magnetic flux, full vector polarimetry of several lines in 
th e present study a nd ma gnetograms of a single line in the case 
of iBerger & Title! yOOl), which result in different noise levels 
and polarimetric sensitivities, and/or (2) the lower spatial reso- 
lution of our observations. 



4.8. Dependence of the temperature on magnetic flux and 
field strength 

A slight reduction in the temperature of the magnetic interior 
with increasing magnetic flux is found, as can be seen in the 
lower left panel of Fig. [18] A similar correlation exists between 
the temperature of the magnetic component and the field strength 
(lower right panel of Fig. [T8l . The temperature of the field-free 
surroundings does not show any correlation with either the mag- 
netic flux or field strength. 



4.7. Dependence of G-band contrast on temperature 



The temperature stratification of the BPs is derived from the in- 
version, both for the magnetic interior and the immediate field- 
free surroundings. To facilitate comparisons, we consider the 
mean temperature of each atmospheric component in the optical 
depth range from log t = to log t = -2. The relation between 
the G-band contrast and the temperature of BPs is displayed in 
the upper left panel of Fig. [18] The contrast shows no clear corre- 
lation with the temperature of either the magnetic interior or the 
field-free surroundings. However, the temperatures of the two 
components show systematic differences: the histograms dis- 
played in the upper right panel reveal that the magnetic atmo- 
sphere is hotter than the immediate field-free surroundings by 
about 1000 K at equal optical depth. The mean temperature is 
5800 K for the magnetic component and around 4800 K for the 
field-free surroundings. 



4.9. The size of G-band BPs 

The quantity / mag = (1 - B)f represents the fraction of a single 
VTT pixel occupied by magnetic fields, as estimated by the in- 
version code. It can be converted to the effective diameter D of a 
circular structure with the same area. The upper limit of D for a 
single pixel is 292 km; however, some BPs cover n VTT pixels. 
The size of each BP in the VTT maps is thus calculated from 
«(/ mag )Bp. In this derivation it is implicitly assumed that only a 
single source is generating the polarized light, even if the mag- 
netic filling fraction only gives the area inside the pixel occupied 
by all fields. For the VTT, no lower size limit exists, because the 
magnetic filling factor can go smoothly towards zero. For the 
DOT images, the number of bright pixels belonging to a single 
BP structure is converted to an effective diameter. The minimum 
effective diameter for any BP in the DOT filtergrams is then 60 
km (= 1 DOT pixel). 

Figure[l9]displays the histogram of BP diameters as inferred 
from both data sets under these assumptions. The maximum di- 
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ameter is 500 km in the VTT maps and 700 km in the DOT 
filtergrams. We find a distribution with a mean value of 150 km 
for the VTT and 210 km for the DOT, with few BPs with diam- 
eters above 300 km. However, the mean value may be mislead- 
ing due to the skewness of the distribution. For comparison, the 
overplotted log-normal distribution would indicate D ~ 90 km 
as the most probable BP diameter. In the DOT filtergrams, there 
is an increase of the relative frequency of small BP structures 
extending over 1 or 2 pixels (60-100 km). These structures are 
very likely artifacts introduced by the data processing and the 
alignment procedure of DOT and VTT data, where the DOT im- 
age closest in time is used without a smooth transition from one 
DOT image to the next. In principle, these smallest BPs should 
have been rejected as unphysical, because the DOT resolution 
of Of.' 2 forces the minimum diameter to be around 170 km. Our 
manual removal of short contour lines thus seems to have been 
insufficient, but fortunately only a small fraction of the BPs is af- 
fected. Note that a more sophisticated algorithm for the identifi- 
cation of BPs would also most probably have resolved the larger 
structures into chains of smaller BPs. The BP sizes found here 
are comparable to the typical wid t h of in terg ranular lanes and th e 
BP sizes given by lBerger & Title! d200lh and lWiehr et alJ(l2004h . 
respectively. 

The lower panel of Fig. [19] compares the effective diame- 
ters of BPs inferred from the VTT and the DOT observations. 
The BP sizes derived from the former are systematically smaller, 
amounting to only 78 % of the DOT diameter on average. The 
smaller size may be related to the spatial resolution of the VTT 
spectra. The magnetic filling fraction is underestimated, as part 
of the polarized light is scattered out of the VTT pixel by seeing. 



5. Summary and discussion 

The analysis of cospatial and simultaneous G-band images and 
spectropolarimetric measurements reveals different aspects of 
both the statistics of BP fields and the relation between G-band 
intensities and magnetic properties of the flux concentrations. 
We find the following properties from a statistical analysis of 
447 identified BPs: 

- 94% of the BPs are cospatial with polarization signals above 
the noise level. 

- The magnetic field strengths of BPs range from 500 to 
1400 G. The field strength distribution is rather flat within 
this range. 

- The total magnetic flux of BPs ranges from to up to 4 x 
10 17 Mx, with an average flux density per pixel of around 50 
G. 

- The distribution of magnetic field inclinations to the LOS, 
hence to the solar surface, indicates that not all BP fields are 
vertical. 

- The BPs observed in a sunspot moat exhibit redshifts in 
the magnetized plasma. The Stokes V zero-crossing veloc- 
ity derived from the IR lines is around 500 m s _1 , while that 
derived from the visible lines is approximately 400 m s -1 . 
The magnetic component of the inversion shows a redshift 
of 260 m s on average. Moreover, the profiles of visible 
and IR lines exhibit amplitude and area asymmetries that in- 
crease with the magnetic velocity. 

- The magnetic interior of BPs is about 1000 K hotter than the 
field-free surroundings at equal optical depth. The tempera- 
ture of the magnetic component is slightly reduced for large 
magnetic fluxes. 
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Fig. 19. Top: Histogram of BP effective diameters from the VTT 
(thin black line) and the DOT (thick line). Overplotted is a log- 
normal distribution. Bottom: Effective diameters from the DOT 
vs those from the VTT. The dashed line has unity slope and 
60 km offset. 



- The area covered by individual BPs in the DOT filtergrams 
is roughly consistent with the size of BPs inferred from the 
inversion of the visible and infrared lines. The average ef- 
fective diameter of BPs is around 100 to 150 km, with few 
structures larger than 300 km. 

For the relations between G-band intensity and magnetic 
properties of the flux concentrations, we find that: 

- The G-band contrast of BPs increases slightly with the mag- 
netic field strength. 

- The contrast only slightly depends on the magnetic flux per 
pixel. It increases with flux for values below 5 x 10 16 Mx and 
decreases with flux for values above 10 17 Mx. 

- The G-band contrast of BPs scales with the inclination: the 
smaller the angle between the magnetic field and the line of 
sight, the higher the G-band contrast. 

Most of these observational findings agree with the expected 
properties of small-scale magnetic flux elements. First, magnetic 
elements in horizontal temperature equilibrium have lower gas 
densities, in order to maintain lateral pressure balance with the 
external field-free medium. Due to the reduced gas density, a 
given optical depth corresponds to lower geometrical heights 
in the magnetic interior than in the field-free surroundings: in- 
side the magnetic flux concentration, deeper and hence hot- 
ter layers are seen. The downward shift of the optical depth 
scale depends on the mag netic field strength, which enters the 
pressure balance equation dSpruitHl976HSteiner & Stenfloll l990: 
iBellot Rubio et al l l2000t ISchussler et alj 120031: IShelvag et al.l 
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Fig. 20. Temperature of the magnetic atmosphere vs its filling 
fraction. For large filling factors, a reduction of the temperature 
is seen. 

120041) . Second, a G-band BP essentially outlines a deficit 
in the abundance of the CH molecule. Under thermody- 
namic equilibrium conditions, the amount of CH in the so- 
lar atmosph ere depends both on the gas density and on the 
temperature (|Sanchez Almeida et al. 2001b; ISteiner et aTl feOOl ; 
Langhans et al.ll2.Q02b . A strong enhancement of the G-band in- 
tensity contrast with respect to nearby continuum wavelengths 
indicates a weakening of the CH lines due to increased thermal 
dissociation. 

We want to explicitly point out the relation between some 
of our observational findings and the modeling of flux concen- 
trations. The temperature difference between the magnetic and 
field-free components of the inversion at the same optical depth 
(Fig. 18) agrees with a shift of the optical depth scale in the 
presence of magnetic fields. The G-band contrast is found to in- 
crease by around 0.2 for fields from 0.4 to 1.5 kG, suggesting a 
stronger shift of the optical depth scale for higher field strengths. 
How ever, the trend falls s hort of the predictions from simulations 
(e.g. lShelvag et al.1 l2004, their Fig. 4), which suggest a contrast 
increase greater than 0.5 over the same range. 

The slight decrease in the G-band contrast for magnetic 
fluxes larger than 10 17 Mx is at first not predicted by models of 
small flux elements. However, we suggest that it is actually re- 
lated to a temperature effect. Since there is an upper limit to the 
field strength on the order of 1.5 kG (cf. Fig. [8}, the amount of 
flux must be proportional to the diameter of the flux concentra- 
tion. The volume inside the flux element to be heated increases 
faster with the diameter of the flux concentration (oc r 2 dh) than 
does the interface with its surroundings (oc rdh). Thus, larger 
magnetic flux should le ad to larger areas, hence to reduced tem- 
peratures (ISpruitlll976l) . While this argument is rather indirect, 
the inversion results allow us to directly check the relation be- 
tween temperature and area. Figure [20] displays a scatter plot of 
the temperature of the magnetic component of the inversion vs 
the magnetic filling fraction, which is proportional to the area. 
In this plot, we see a clear trend toward a reduction of the tem- 
perature with filling factor, hence with area. 

Perhaps the most surprising result of our analysis is the de- 
pendence of the G-band contrast on the magnetic field inclina- 
tion to the LOS displayed in the middle panels of Fig. [TT] We 
note that the dependence seems to have the form 

b 

C = a + , (11) 

siny 

where y represents the field inclination and a and b are constant 
parameters. A very simple model of a vertical flux concentration 
embedded in field-free surroundings (Fig.l2"TTi shows that the sec- 



Fig. 21. Sketch of an isolated flux concentration embedded in 
field-free surroundings. The magnetic field is along the vertical 
direction and the LOS is inclined by an angle y with respect to 
it. The path length of a ray inside the flux concentration, AR, 
depends on both y and the thickness d of the tube. 

ond term on the right hand side of Eq. (fTTT > is proportional to AR, 
the distance inside the flux concentration traveled by a ray that 
hits the magnetic element at an angle y. In fact, AR = d/ sin y. 
The absorption along the LOS would then be oc AR ■ p(B), with 
the density p(B) inside the flux concentration being lower than 
the density outside it. Thus, AR would finally translate into a 
downward shift of the optical depth scale for lines of sight that 
pass through the magnetic element, compared with others that 
miss it. Note that AR is limited by a AR max , the depth where the 
atmosphere gets opaque even for y — deg. We have overplotted 
three curves for different values of a and b in the middle panels 
of Fig. [TT] where the middle curve (b - 200) at least roughly 
agrees with the observed inclinations and L/V-values. 

6. Conclusions 

We have derived the thermal, magnetic, and kinematic proper- 
ties of G-band bright points in the moat of a regular sunspot 
from infrared and visible spectral lines observed at the German 
Vacuum Tower Telescope on Tenerife. The BPs were identified 
in cospatial diffraction-limited filtergrams taken with the Dutch 
Open Telescope on La Palma. The Stokes profiles of the infrared 
and visible spectral lines were inverted simultaneously using 
a two-component model atmosphere with height-independent 
magnetic fields and line-of-sight (LOS) velocities. In addition, 
line parameters were extracted from the observed Stokes profiles 
for an estimate of field strengths, field inclinations, and LOS ve- 
locities. 

We conclude that G-band brightenings are caused by con- 
centrated magnetic fields in more than 90% of the cases. The 
minimum requirement seems to be a field strength of at least 
500 G. However, the G-band BPs show a variety of magnetic 
and kinematic properties, as suggested by the broad range of 
values we find in all physical quantities (B, y, O, etc). Several 
patches of magnetic flux can be traced in the polarimetric data 
during the full 1-hour time series, which show up as BP only 
part of the time. Together with the broad range of field strengths, 
we think this indicates that not all BPs are cospatial with sta- 
ble kG flux tubes. We find a stronger dependence of the G-band 
brightness on a geometrical effect, namely on the inclination of 
the magnetic field to the LOS, than on most of the other phys- 
ical quantities. We conclude that for an accurate description of 
the generation of G-band BPs or, more generally, the intensity in 
the G band, it is necessary to develop models of flux concentra- 
tions and their surroundings in at least two dimensions, to take 
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into account the geometry of the field, its strength, and the actual 
viewing angle of the observations. 

Even if our data corresponds to the moat of a sunspot, we 
believe that the magnetic properties of the observed BPs are not 
strongly influenced by the presence of the spot. The distance to 
the spot boundary is several arcseconds in most cases. Thus, the 
results of this paper might also apply to BPs observed in other 
active and non-active regions. 

Finally, we stress that the use of G-band brightness enhance- 
ments as proxies of magnetic fields may miss part of the solar 
magnetic flux, as fields with strengths below 500 G remain un- 
detected. 
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Appendix A: Spatial resolution 

For an estimate of the spatial resolution, we calculated the 
Fourier transform of the intensity maps of the first scan across 
the sunspot. We then integrated the power of the Fourier trans- 
form over rings corresponding to increasing spatial frequencies. 
The inverse of the spatial frequency corresponds to a spatial 
scale. The resolution limit was estimated as the point, where 
the power spectrum levels off into the constant noise contribu- 
tion. For the G-band data from the DOT, we find a resolution of 
around 0/25 for the original DOT map with O/071 square pixels 
(cf. Fig. IA. lb . When the data is resized to the VTT resolution 
of 0737 square pixels, the power is still conserved. For the po- 
larimetric data in either infrared or visible continuum, we find a 
resolution of 170 with this method. 



Appendix B: Alignment of data sets 

The alignment procedure uses the TIP pixel as a reference sys- 
tem: for each position (x,y)xip in the TIP map, where x corre- 
sponds to the scan step and y to the position along the slit, the 
algorithm looks for cospatial points in the visible and UV chan- 
nel of POLIS and the DOT filtergrams, respectively. The TIP 
data have a spatial sampling of 0735 per pixel. The data from 
the visible channel of POLIS have a finer sampling of 07145, 
which is degraded to the TIP resolution by the linear interpo- 
lation scheme described in the next section. The data from the 
VTT show displacements due to the actual position of the TIP 
and POLIS cameras and an additional var iable offset caused by 
differential refraction (e.g. lReardonll2006l) in the Earth's atmo- 
sphere. 

B. 1 . Alignment between IR and visible/UV data 

Calculation of shifts. To find cospatial points in the TIP and 
POLIS data, the continuum intensity maps of each repetition of 
the scan are correlated separately. For the UV channel, a map 
of the line wing intensity of Ca n H is used, as it shows photo- 
spheric structure (cf. Fig.[TJ. The resulting shift values along and 



Fig. B.l. Displacements between the visible and infrared inten- 
sity maps (+), and between the UV and infrared intensity maps 
(❖), in arcsec. For each repetition of the scan across the spot dis- 
played in Fig.[T] the shifts perpendicular to the slit (Ax; top) and 
along the slit {Ay; bottom) have been calculated through cross 
correlation. The dotted and dashed lines give the shift values 
used for the alignment of each scan step. 

perpendicular to the slit for the different times of the eight repe- 
titions allow a curve to be fitted to the displacements (Fig. IB. It . 
The best-fit curve is used to determine the shifts in x and y 
needed to align the TIP and POLIS observations for each scan 
step (taken every 6 seconds). 

Application of shifts. The spatial sampling of 0729 in the UV 
channel of POLIS is similar to that of TIP. For this reason, the 
shift (Axuv, Ayt/y) is simply added to the position (x,y)jip, and 
the closest pixel is taken: 

(x, y) uv = round [(x, y) T ip + (Ax uv , Ay uv )] , (B.l) 

where the rounding automatically selects the nearest cospatial 
pixel. 

The spatial sampling in the visible channel of POLIS is al- 
most three times better than that of TIP. The cospatial Stokes 
profiles are then retrieved as a weighted average of the visible 
profiles along and perpendicular to the slit. First, the cospatial 
position (x, y) v ; s is given by 

(x,y) vis = (x,y) TlP + (floor[Ax vis ], Ay vis ), (B.2) 

where floor[Ax v ; s ] returns the largest integer smaller than x v ; s . 

The TIP pixel partly covers several POLIS pixels from two 
different scan steps, x v ; s and x v j s + 1 (Fig. |B.2t . The cospatial 
Stokes profile p(A, x, j) v ; s is calculated as the weighted average 

p(A,x,y) vis = + — , (B.3) 

Li «« Li hi 

where i = 1, . . .,3 or i = 1, ... ,4. The weights a, and b t are set 
according to the fraction of the POLIS pixels covered by the TIP 
pixel and verify the condition + bf) = 1. 



18 



C. Beck et al,: Magnetic properties of G-band bright points 



x+1 



scanning direction 



0.35" 







p 4 


















p 3 




(x, 








(x, 






p 2 


















Pi 

















| slit 
^0.145" 



0.35" 



r, 



0.35" 



Fig. B.2. Calculation of the visible profile that is cospatial to the 
infrared Stokes profile at (x,v)tip- Two cases are possible: the 
TIP pixel covers three (left) or four (right) POLIS pixels along 
the slit. The profiles p, from the two scan steps x v \ s and jt v j s + 1 
are averaged with weights corresponding to the fraction of the 
POLIS resolution elements covered. 

B.2. Alignment between DOT and VTT data 

In this case, cospatial and cotemporal positions have to be found. 
A cotemporal map from the DOT time series ("pseudo-scan") is 
first constructed and then spatially aligned with the VTT data. 

B.2.1. Cotemporal DOT map 

The DOT pseudo-scan is the image that would result from step- 
ping the POLIS/TIP slit across the DOT FOV. To create this map, 
the POLIS slit-jaw images are used. They display a 100"x 100" 
FOV of the solar surface centered on the slit. Each slit-jaw image 
is correlated with the DOT blue continuum filtergram closest in 
time, degraded to the same spatial resolution. From the correla- 
tion one finds the shifts in x and y required to align the images 
from the two telescopes. The POLIS slit visible in the slit-jaw 
image is superimposed on the coaligned DOT filtergram. A slice 
of 5 pixel width (= O^SS) is taken from the DOT images at 
the slit position and placed accordingly in the pseudo-scan map. 
Examples of the resulting pseudo-scan maps for the three DOT 
channels are displayed in Fig.Q] 

B.2.2. Cospatial DOT map 

The final alignment of the DOT images is achieved by taking a 
25x25 pixel subfield of the TIP map around each pixel (x,y)xip, 
the corresponding area from the DOT pseudo-scan continuum 
map, and correlating the two subfields. This procedure yields 
the position of a 5x5 pixel area in the DOT maps that is cospa- 
tial with the TIP pixel. The area is cut out from the three DOT 
pseudo-scan maps and placed in new maps accordingly. Fig. [3] 
displays the coaligned maps for the first repetition of the scan. 
The flow chart displayed in Fi g . IB . 3 1 summarizes all the steps we 
follow to align the various data sets. 



Co-spatial 
information 



vis. profile 




interpolation 
/ shift vis 



(x,y) in IR map 



Reference 
position 



subfield x+dx 

y+dy 



correlation 



subfield x+dx 

y+dy 



DOT 'pseudo 
scan 1 from 
slices along 
POLIS slit 



aligned DOT data 



Fig. B.3. Alignment procedure. For each position (x,y)xip in the 
TIP data set, cospatial POLIS profiles are found by applying ap- 
propriate shifts to the visible and UV channels. To align the DOT 
filtergrams, the POLIS slit-jaw image and the DOT continuum 
image closest in time are correlated and used for the creation of 
a pseudo-scan map. The final DOT/VTT alignment is done by 
taking a subfield of the TIP map around (x, y)np, the correspond- 
ing subfield from the DOT pseudo-scan, and then correlating the 
subfields. 



A single polarization patch usually contains a few BPs, and the 
majority of polarization patches show BPs at least in one of the 
repetitions. Most magnetic signals remain visible during the full 
observation run (1 hour). This agrees with the interpretation that 
brightenings in the G band are due to the presence of magnetic 
fields, while the high temporal variability of the BPs is due to 
changes in the magnetic field configuration. 

The sunspot canopy (visible near the upper edge of the FOV) 
clearly extends beyond the white-light boundary of the spot. The 
canopy is not an area of enhanced G-band intensity because its 
magnetic fields are very inclined and do not lead to a sufficient 
downward shift of the optical depth scale. 



Appendix C: Coaligned maps 

Coaligned DOT and VTT maps are displayed in Figs. IC.ll and 
IC.2l for the eight repetitions of the scan. We show the aligned 
DOT G-band, the total integrated polarization, the IR continuum 
intensity, and the DOT BP mask. For the first repetition of the 
scan, the VTT mask of identified BPs is displayed instead of 
the DOT BP mask. The area used to extract the NBP sample is 
also indicated for the first repetition. A visual inspection shows 
that almost all features identified as G-band bright points (red 
contours) are cospatial with polarization signals above the noise. 
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Fig. C.l. Coaligned maps from the VTT and the DOT for the first four repetitions of the scan. Clockwise, starting from top left of 
each 4x4 subpanel: G-band intensity, total integrated polarization, DOT BP mask, and IR continuum intensity. Blue contours outline 
enhanced polarization signals and red contours the selected BP areas. In the DOT BP mask, the BPs are color-coded according to 
their number. The lower right map of the first 4x4 subpanel is the VTT BP mask, not the DOT BP mask. In this map, the uniform 
pink area marks the spatial points used to extract the NBP sample. They are located outside the canopy, and at least 3 pixels (1'.'05) 
away from any BP. 



Fig. C.2. Same as Fig. IC. II for the last four repetitions of the scan. 
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